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Cycle  of  vertical  and  horizontal  mixing  in  a  shallow  tidal  creek 

Paul  McKay1,2  and  Daniela  Di  Iorio1 
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[i]  An  experimental  study  of  vertical  mixing  and  along-channel  dispersion  parameterized 
in  terms  of  horizontal  mixing  near  the  mouth  of  the  Duplin  River  (a  tidal  creek  bordered 
by  extensive  intertidal  salt  marshes  on  Sapelo  Island,  Georgia)  was  carried  out  over 
several  spring/neap  cycles  in  the  fall  of  2005.  Vertical  mixing  is  modulated  on  both  M4 
and  fortnightly  frequencies  with  maximum  turbulent  stresses  being  generated  near  the 
bed  on  periods  of  maximum  flood  and  ebb  and  propagating  into  the  water  column 
showing  a  linear  dependence  with  depth.  Values  are  significantly  greater  on  spring  tide 
than  on  neap.  Horizontal  mixing  evaluated  by  salt  fluxes  is  driven  and  dominated  by 
tidal  dispersion,  which  is  also  modulated  by  the  fortnightly  spring/neap  cycle.  Net  export 
of  salt  from  the  lower  Duplin  is  shown  to  be  due  to  residual  advection  modified  by 
upstream  tidal  pumping.  The  tidal  dispersion  coefficient  exhibits  a  pulsating  character 
with  greater  values  on  spring  tide  followed  by  smaller  values  on  neap  tide. 

Citation:  MeKay,  P.,  and  D  Di  Iorio  (2010),  Cycle  of  vertical  and  horizontal  mixing  in  a  shallow  tidal  creek,  J.  Geophys.  Res.,  115 , 
CO  1004,  doi:  1 0. 1 029/2 008JC005204. 


1.  Introduction 

[2]  Efforts  to  understand  the  nature  of  the  transfer  of 
materials  and  dissolved  substances  across  the  land-ocean 
boundary,  and  into  the  coastal  ocean,  must  necessarily  focus 
on  the  estuaries  and  tidal  creeks  where  such  fluxes  are 
concentrated.  Since  the  early  work  on  estuarine  dynamics 
[see,  e.g.,  Pritchard ,  1 952,  1 954;  Hansen  and  Rattray ,  1 965] 
vertical  mixing,  largely  dnven  by  bottom  stress,  surface 
stress  and  internally  generated  turbulence,  has  been  known 
to  be  an  important  control  on  the  salt  flux  through  its  action 
on  vertical  stratification.  More  recently  the  spring-neap 
cycle,  with  its  attendant  modulation  of  tidal  energy,  has  come 
to  be  appreciated  as  an  important  control  on  this  vertical 
mixing  itself  [see,  e.g.,  Peters ,  1997;  Chant ,  2002;  Simpson 
et  al. ,  2005;  Chant  et  al. ,  2007], 

[3]  Longitudinal  salt  fluxes  arc  the  result  of  both  advee- 
tive  and  dispersive  processes  [Dyer,  1997],  with  the  advec- 
tive  fluxes  due,  primarily,  to  fresh  water  input  at  the  head 
and  the  dispersive  fluxes  due  to  gravitational  circulation  and 
tidal  processes,  and  thus  maximum  near  the  mouth  of  an 
estuary  or  tidal  creek  [Geyer  and  Signed ,  1992].  These 
longitudinal  fluxes  can  vary  greatly  on  seasonal  and  spring/ 
neap  time  scales,  with  changes  in  fresh  water  input  and  tidal 
energy  making  them  difficult  to  predict  [see,  e.g.,  Banas  et 
al ,  2004;  Medeiros  and  Kjerfve ,  2005]. 

[4]  This  paper  presents  the  results  of  a  experiment  con¬ 
ducted  over  several  spring/neap  cycles  near  the  mouth  of  the 
Duplin  River,  a  salt  marsh  tidal  creek  on  the  central  Georgia 

1  Dcpartmcni  of  Mannc  Science,  University  of  Georgia,  Alhcns, 
Georgia,  USA. 

2Now  al  Naval  Research  Laboratory,  Stcnnis  Space  Center,  Mississippi, 
USA. 
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coast.  The  objective  of  this  study  is  to  quantify  the  influence 
of  the  spring/neap  cycle  on  both  vertical  mixing  and 
longitudinal  dispersion,  which  can  be  viewed  as  horizontal 
mixing,  with  the  aim  of  better  understanding  the  relative 
importance  of  the  down  estuary  flux  by  freshwater  flow  and 
the  up-estuary  flux  by  dispersive  processes.  The  mean 
properties  of  the  Duplin  have  been  previously  reported  on 
[see,  e.g.,  Ragotzkie  and  Bryson ,  1955;  Kjerfve,  1973; 
Imberger  et  al.,  1 983]  but  this  is  the  first  study  in  this  region 
to  look  at  both  horizontal  and  vertical  mixing  over  spring/ 
neap  time  scales,  processes  that  are  important  for  controlling 
the  distribution  and  dilution  of  salinity.  As  tidal  mixing  is 
strong  in  this  shallow  estuary  (as  will  be  shown),  it  is 
expected  that  any  stratification  induced  by  straining  pro¬ 
cesses  will  have  a  negligible  effect  and  thus  we  classify 
this  estuary  as  well  mixed. 

[5]  The  study  site  is  a  tidal  creek  located  in  the  marshes 
of  the  Sapelo  Island  National  Estuarine  Research  Reserve 
on  the  central  Georgia  coast  in  the  southeastern  U.S.  (see 
Figure  1).  The  Duplin  is  approximately  13  km  long,  three 
times  the  estimated  tidal  excursion  distance  [Ragotzkie  and 
Bryson,  1955],  and  varies  from  less  than  1  m  wide  at  its 
head,  in  the  intertidal  marsh,  to  200  m  wide  at  its  mouth.  It 
connects  with  Doboy  Sound,  which  in  turn  is  connected  to 
the  coastal  Atlantic  Ocean.  The  upper  reaches  of  the  Duplin 
wind  through  extensive  intertidal  marshes  characterized  by 
mud  flats,  vegetated  primarily  with  Spartina  altemiflora 
and  cut  by  a  network  of  side  creeks  and  channels  of  varying 
sizes,  while  the  lower  reaches  consist  of  a  straight  and  wide 
channel  with  few  side  creeks  bordered  by  more  upland, 
marsh  hammocks  and  creekless  marsh. 

[6]  The  mean  depth  along  the  thalweg  in  the  lower 
Duplin  is  approximately  6.5  m.  The  drainage  area  of  the 
intertidal  marsh  is  approximately  12  km2  [Blanton  et  al., 
2007].  Tidal  flows  are  ebb  dominant  and  this  is  a  common 
characteristic  of  this  class  of  salt  marsh  estuary  [Dronkers, 
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Figure  1.  Chart  showing  the  study  area  on  the  central  Georgia  coast  showing  the  ADCP  mooring  in  the 
lower  Duplin  (asterisk),  the  Marsh  Landing  weather  station  (solid  square),  and  the  GCE- 1 0  site  in  the  upper 
Duplin  (solid  circle)  where  the  GCE-LTER  project  maintains  a  long-term  monitoring  station.  The  offshore 
National  Data  Buoy  41008  (open  square)  provides  the  offshore  wind  forcing. 


1986;  Dyer ,  1997;  Blanton  et  a/.,  2002]  with  a  short  strong 
ebb  flow  and  a  slower  prolonged  flood. 

[7]  The  freshwater  input  in  this  region  is  primarily  from 
submarine  groundwater  input  [Ragotzkie  and  Bryson ,  1 955] 
with  intermittent  surface  runoff  due  to  precipitation.  This 
input,  which  is  ungauged,  establishes  an  along  channel 
salinity  gradient  in  the  lower  reaches  which  can  be  as  much 
as  2.5  per  km  or  more  (as  measured  using  the  praetieal 
salinity  scale)  over  the  length  of  one  tidal  excursion. 
Upstream  of  this  lower  section  the  waters  are  generally  well 
mixed  along  the  channel  and  the  salinity  gradient,  when 
present,  is  much  smaller  and  variable.  There  is  significant 
tidal  salinity  variation  in  the  lower  reaches  of  the  Duplin  as 
water  from  Do  boy  Sound  (influenced  by  the  Altamaha  River) 
is  advcctcd  by  Kjerfve  [1973].  While  there  are  occasional 
periods  of  stratification  in  the  lower  Duplin  (believed  to  be 


related  to  changes  in  Altamaha  River  discharge,  groundwater 
input  and  precipitation),  it  is  generally  vertically  and  laterally 
well  mixed 

[s]  In  this  paper,  section  2  outlines  the  details  of  the 
experiment  and  data  processing  and  section  3  describes  the 
mathematical  framework  for  vertical  and  horizontal  mixing 
and  dispersion  in  estuaries.  The  measured  mixing  and  dis¬ 
persion  results  are  summarized  in  section  4  and  section  5 
discusses  implications  fora  better  physical  understanding  of 
tidal  creek  export  to  the  coastal  ocean. 

2.  Experimental  Program 

[9]  Meteorological  conditions  at  the  offshore  National 
Data  Buoy  41008  station  and  the  Marsh  Landing  weather 
station  adjacent  to  the  deployment  location  (see  Figure  1) 
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Figure  2.  The  coastal  wind  stress  in  terms  of  along-  (positive  is  to  the  north)  and  cross-  (positive  is  to 
the  east)  shelf  components  showing  that  during  this  time  of  year  the  winds  are  predominantly  from  the 
northeast.  The  precipitation  rate  at  the  Marsh  Landing  weather  station  immediately  adjacent  to  the 
mooring  location  arc  also  shown. 


are  shown  in  Figure  2.  Wind  stresses  show  that  the  coastal 
storms  are  predominately  from  the  northeast  which  can 
increase  inundation  in  the  estuary  through  Ekman  transport 
processes  (as  will  be  seen  in  the  residual  sea  level  height  in 


Figure  4)  Sporadic  rain  showers  were  observed  on  year  day 
274  281,  295  and  297,  and  the  region  had  been  recovering 
from  a  drought  in  2001  and  2002.  Some  of  the  decrease  in 
mean  salinity  from  year  day  278  281  (as  will  be  shown)  is 
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Figure  3.  The  depth-dependent  along-ehannel  velocity  («)  is  plotted  together  with  the  surface  and 
bottom  velocities.  The  Reynold’s  stress  (rjp)  is  shown  during  a  spring-neap-spring  period  from  year  day 
277  291.  Measurements,  where  indicated,  are  a  function  of  mab. 
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likely  attributable  to  this  rain  event,  however  this  period  also 
corresponds  with  the  spring  tide,  which  is  believed  to 
increase  fresh  groundwater  flux  into  the  channel  (as  will  be 
discussed). 

[10]  The  field  experiment  in  the  lower  Duplin  River  was 
carried  out  over  a  34  day  period  from  28  September  to  31 
October  2005  and  covered  two  spring  and  two  neap  tides, 
ending  shortly  before  a  third  spring  tide.  Figure  1  shows  a 
hydrographic  chart  of  the  study  area  on  the  Georgia  coast  as 
well  as  the  location  of  our  instruments.  The  instrument 
package  was  deployed  during  the  September  2005  Georgia 
Coastal  Ecosystems  (GCE)-Long  Term  Ecological  Research 
(LTER)  survey  cruise  with  the  R/V  Savannah.  A  4  beam 
1200  kHz  broadband  acoustic  Doppler  current  profiler 
(ADCP)  was  deployed  on  a  heavy  custom  built  pyramidal 
bottom  mount  mooring  frame  in  the  center  of  the  channel  in 
water  with  a  mean  depth  of  approximately  6.5  m.  The  ADCP 
operated  in  burst  mode,  sampling  current  velocity  profiles  at 

2  Hz  for  5  minutes  every  30  minutes,  and  logged  every  ping 
in  beam  coordinate  mode  with  0.25  m  depth  bins.  The  first 
sample  was  centered  1 .0  m  above  bottom  when  taking  into 
account  blanking  distance.  The  I  2  beam  pair  was  aligned 
along  the  main  axis  of  the  channel 

[n]  For  each  beam,  the  velocity  data  at  the  top  bin  and 
bottom  two  bins  were  deleted,  to  remove  both  surface  noise 
and  ringing  near  the  transducer.  Each  beam  was  then  further 
processed  by  removing  all  ping  data  with  a  percent  good 
(PG)  of  less  than  90%  as  reported  by  the  internal  diagnostics 
of  the  ADCP.  For  an  ADCP  logging  every  ping  this  is  a 
measure  of  the  echo  strength  of  the  returned  signal.  This 
resulted  in  the  removal  of  one  complete  five  minute  ensemble 
on  year  day  289  as  well  as  the  removal  of  scattered  other 
values,  generally  near  the  surface,  which  constituted  less 
than  1%  of  the  data.  With  the  exception  of  the  removed 
ensemble,  deleted  beam  velocities  were  replaced  with  a 
linear  interpolation  between  the  values  in  the  bins  above 
and  below.  Finally,  for  each  depth  bin,  spikes  in  the  data  were 
rejected  that  were  three  or  more  standard  deviations  away 
from  the  ensemble  average  and  were  replaced  with  an  inter¬ 
polated  value.  This  also  removed  a  small  amount  of  data. 

[12]  The  along  channel  (w)  and  cross  channel  (v)  velocities 
are  calculated  from  the  individual  ADCP  pings  as 

b2-b1  b4-b> 

2  sin  O'  V  2  sin  0 

where  B,  is  the  along  beam  velocity  for  each  beam  (/  =  1, 
2,  3,  4),  the  I  2  beam  pair  is  aligned  with  the  channel,  the 

3  4  beam  pair  is  aligned  cross  channel  and  0  =  20°  is  the 
beam  angle  (for  the  RD1  ADCP  used)  [Di  lorio  and  Gargett , 
2005].  The  instantaneous  u  and  v  velocities  were  then 
averaged  for  each  five  minute  ensemble  to  get  the  mean  flow. 
To  resolve  the  surface  and  bottom  flows  each  velocity  profile 
was  extrapolated  to  the  surface,  over  a  distance  of  0.5  m,  by 
assuming  zero  shear  at  the  surface,  and  extrapolated  to  the 
bottom  assuming  a  log  layer  profile  over  a  1 .5  m  distance. 

[13]  Compass  and  pitch/roll  sensors  show  that  during  this 
deployment  the  instrument’s  heading  changed  0.5°  over  34 
days  while  pitch  and  roll  started  out  at  less  than  3  and  1.5°, 
respectively,  from  the  vertical  with  both  settling  to  just  over 
1 0  over  the  first  1 0  days  of  the  deployment  and  remaining 
constant  for  the  rest  of  the  deployment.  While  stress 


estimates  are  very  sensitive  to  instrument  alignment  with 
the  vertical,  Lu  and  Lueck  [  1 999]  showed  that  errors  in  the 
stress  estimate  will  be  small  for  such  small  deviations.  Our 
coordinate  system  adheres  to  the  estuarine  convention  such 
that  x  is  the  along  channel  direction  and  is  positive  out  of 
the  creek. 

[m]  Microeat  eonductivity-temperature-depth  (CTD) 
meters  were  deployed  with  the  ADCP  mooring  with  a 
bottom  sensor  mounted  to  the  ADCP  frame  approximately 
0.5  meters  above  bottom  and  a  surface  sensor  attached  to  a 
line  tethered  to  a  surface  marker  float  and  deployed  0.5  m 
below  the  surface  In  addition  the  GCE-LTER  project 
maintains  a  long-term  CTD  station  in  the  upper  Duplin  at 
mid  depth.  All  CTDs  sampled  at  15  minute  intervals.  The 
bottom  Microcat  salinity  data  became  unreliable,  due  to 
biofouling  of  the  sensor,  on  the  fifth  day  of  the  deployment. 
Prior  to  this  instrument  failure,  the  surface  and  bottom 
measurements  indicated  that  the  system  was  vertically  well 
mixed  for  both  temperature  and  salinity  at  all  times  with  the 
differences  showing  a  random  distribution  and  being  below 
the  noise  threshold  of  each  sensor.  As  this  agrees  with  the 
authors’  observations  during  two  previous  mooring  experi¬ 
ments  at  this  location  in  the  Duplin  River  [McKay  and 
Di  lorio ,  2008]  we  have  chosen  to  treat  the  surface  values 
as  being  indicative  of  the  entire  water  column.  While  rain 
events  likely  introduced  some  short-lived  stratification  to 
the  water  column  we  were  not  able  to  measure  this.  Salinity 
will  be  reported  as  a  unitless  quantity  in  accordance  to  the 
practical  salinity  scale  adopted  in  1978. 

3.  Methods 
3.1.  Vertical  Mixing 

[15]  In  a  sheared  flow,  vertical  mixing  results  from 
bottom  boundary-induced  turbulence  or  breaking  internal 
waves  (Kelvin-Helmholtz  instabilities)  depending  on  the 
strength  of  the  stratification.  For  the  weakly  stratified  ease. 
Jay  and  Smith  [1990]  show  that  when  the  shear  is  large 
enough  relative  to  the  stratification  then  bottom  boundary- 
induced  turbulent  mixing  extends  throughout  most  of  the 
water  column  and  causes  vertical  exchange  of  properties. 
Any  tidal  asymmetries  would  then  come  from  the  mean 
flow  as  the  velocities  on  ebb  are  greater  because  the  shear 
and  flow  speed  are  augmented  by  the  mean  river  flow.  The 
barotropie  tidal  adveetion  of  the  density  field  dominates  in 
weakly  stratified  estuaries,  as  the  gravitational  restoring 
force  is  weak  and  therefore  baroclinie  motions  do  not  exist. 

[16]  It  is  common  to  treat  turbulent  mixing  (vertical  flux 
of  horizontal  momentum  and  salt)  as  a  gradient  process 
parameterized  against  the  mean  velocity  shear  and  salinity 
gradient  with  a  vertical  eddy  viscosity,  A2 ,  and  vertical 
scalar  diffusivity,  Az,  respectively,  (with  units  of  m2  s  l) 
such  that 


(2) 

-  Os 

Oz 

(3) 

The  along  channel  and  vertical  velocities  are  u  and  >v, 
respectively,  and  the  subscript  z  is  the  vertical;  the  ovcrlinc 
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represents  a  time  averaged  quantity  and  primed  values 
represent  the  turbulent  fluctuations  from  the  time  average 
l Dyer ,  1997]. 

[17]  Since  the  work  of Lohrman  et  ai  [1990],  as  elaborated 
by  Stacey  et  ai  [1999a,  1999b]  and  Lu  and  Lueck  [1999],  it 
has  become  common  to  estimate  vertical  eddy  viscosity  (Az), 
and  hence  momentum  diflusivity,  using  a  four  beam,  high- 
frequency  broadband  acoustic  Doppler  current  profiler 
(ADCP).  This  method,  known  as  the  variance  method,  £l]ows 
the  direct  estimation  of  the  Reynolds  stress,  rjp  -  — u'v 1/,  by 
comparing  the  velocity  variances  of  opposing  pairs  of 
beams.  This  stress  is  determined  using  the  relation 


—iSwJ  = 


Bh  ~  Bl 


1/ 


2  sin  20  ’ 


(4) 


where  B\f  =  B]  —  B\ ,  B\  is  the  velocity  along  the  beam  one 
direction  and  the  subscript  /  implies  fluctuations  [Di  lorio 
and  Gargett,  2005].  B\^  is  similarly  defined  but  opposite 
from  beam  one.  The  vertical  viscosity  Az is  then  the  ratio 
between  the  time  varying  stress  and  the  observed  vertical 
shear. 

[18]  In  energetic  tidal  channels  having  vertically  well- 
mixed  conditions  and  hence  a  gradient  Richardson  number 
tending  to  0,  the  turbulent  Prandtl  number  AJKZ  =  0.7 
[Stacey  et  ai ,  1999a].  In  the  presence  of  stratification, 
vertical  mixing  has  to  act  against  a  density  gradient  and 
both  the  eddy  viscosity  and  diflusivity  are  reduced.  How¬ 
ever,  as  momentum  is  transferred  more  readily  than  mass  in 
such  a  case,  Kz  is  attenuated  more  than  is  Az  and  the 
turbulent  Prandtl  number  increases. 

[19]  The  variance  method  requires  the  flow  to  be  hori¬ 
zontally  homogeneous  such  that  there  is  no  variation  in  the 
turbulent  statistics  over  the  separation  distance  between 
corresponding  bins  in  each  beam  pair  [Lu  and  Lueck , 
1999].  This  condition  is  met  by  our  deployment  location 
being  in  the  middle  of  the  channel  in  a  region  of  nearly 
constant  cross  section.  The  sampling  rate  of  2  Hz  and 
vertical  bin  size  of  0.25  m  allows  us  to  resolve  the  smaller 
eddies  in  the  flow  which  are  involved  in  vertical  momentum 
transfer  [Rippeth  et  ai ,  2002].  The  5  minute  ensemble  time 
represents  a  compromise  between  concerns  of  instrument 
battery  life,  the  need  for  a  statistically  significant  sample  size 
and  the  need  for  quasi -station ary  conditions  during  the  sample 
period.  With  a  2  Hz  sampling  rate,  a  5  minute  sampling  period 
(600  samples)  provided  a  statistically  significant  ensemble 
average  of  the  velocity  as  the  estimated  standard  error  was 
low  (less  than  1  cm/s).  In  a  flow  dominated  by  semidiurnal 
tides,  the  5  min  period  was  also  short  enough  to  insure  quasi¬ 
stationary  conditions  during  the  entire  ensemble. 

[20]  Stress  profiles  were  calculated  for  each  ping  in  the 
record  and  then  averaged  for  each  five  minute  ensemble  to 
give  one  stress  profile  for  each  half  hour  period.  The 
statistical  reliability  of  the  Reynolds  stress  estimates 
increases  with  the  square  root  of  the  number  of  samples  per 
ensemble  according  to  the  relation  given  by  Williams  and 
Simpson,  [2004] 


a 


2  _ 

R  — 


^-4)2 

N  sin2  20 


(5) 


where  o\  is  the  mean  squared  variability  of  the  Reynolds 
stress  estimate,  7  is  a  factor  depending  on  the  covariance  of 
the  individual  velocity  values,  <jn  is  the  instrument  noise 
level,  Bz  the  mean  square  value  of  the  turbulent  fluctuations 
and  N  is  the  number  of  samples  per  ensemble.  From 
equation  (5)  and  following  the  method  of  Williams  and 
Simpson  [2004],  we  can  estimate  a  threshold  value  for 
detectable  stress  by  examining  low  flow,  when  B,/goes  to 
zero  and  7  to  one,  and  taking  a  value  for  instrument  noise  of 
on  =  0.017  m  s-1  from  the  instrument  manufacturer.  This 
then  gives  an  approximate  theoretical  minimum  measurable 
stress  rfp  =  2  x  10  s  m2  s  2. 

[21]  To  estimate  the  vertical  eddy  viscosity,  Az,  it  is 
necessary  to  compute  the  vertical  shear  in  the  along  channel 
velocity,  dii/dz ,  in  the  region  where  stresses  are  resolved  by 
the  ADCP.  As  small  variations  in  u  between  depth  bins  can 
cause  large  swings  in  the  value  of  duldz ,  as  calculated  using 
numerical  differentiation  techniques,  each  velocity  profile 
was  smoothed  by  fitting  with  a  polynomial  and  the  gradient 
taken  by  differentiating  that  polynomial.  As  the  shape  of  the 
velocity  profile  changes  through  the  tidal  cycle  from  being 
nearly  linear  to  showing  slight  curvature,  each  profile  was 
fit  with  both  first  and  second-order  polynomials;  the  fit  with 
the  highest  r2  value  was  selected  for  the  differentiation.  The 
velocity  profile  was  generally  linear  with  a  maximum  at  the 
surface. 

[22]  The  vertical  mixing  time  during  the  tidal  cycle  is 
estimated  from  a  tidal  and  depth  average  of  the  vertical 
mixing  quantity  Kz  [Lewis,  1997],  as 

//2  //2 

*  —  \  %  /-  \  ’  W 

8  {K:)  8  (Az) 

where  t  is  the  approximate  time  for  complete  vertical  mixing 
(in  seconds),  //  the  tidal  mean  water  depth,  and  ()  denotes 
tidal  averaging  and? denotes  depth  averaging.  A  turbulent 
Prandtl  scaling  of  1  (Kz  ~  Az)  is  assumed  and  that  it  is 
considered  stationary  within  each  tidal  period.  This  estimated 
mixing  time  is  based  on  the  estimated  time  required  for  a 
concentrated  substance  in  the  middle  of  the  water  column  to 
assume  a  Gaussian  concentration  distribution  with  depth 

3.2.  Decomposition  of  Salt  Fluxes 

[23]  Along  channel  horizontal  mixing  in  an  estuary  is  a 
result  of  the  sum  of  several  processes  on  widely  distributed 
spatial  and  temporal  scales:  from  turbulent  to  tidal  and 
fortnightly  scales  The  interaction  between  estuarine  tides, 
bathymetry  and  irregularities  at  the  channel  edge  can  result 
in  dispersive  fluxes  that  are  due  to  tidal  pumping/trapping, 
vertical  and  lateral  shear  dispersion,  and  the  estuarine 
circulation  [Fischer  et  ai,  1979].  Of  these  mechanisms  tidal 
dispersion,  the  generally  upstream  net  transport  due  to 
correlations  between  tidal  depth,  velocity  and  salinity,  is 
often  dominant  in  well  mixed  estuaries  [ Kjerjve ,  1986].  The 
presence  of  vertical  shear  and  stratification  [Bowden,  1965] 
and  a  generally  much  smaller  effect  of  lateral  shear  and 
salinity  variability  [Rattray  and  Dworski ,  1980]  can  act  to 
modulate  this  tidal  pumping  mechanism,  especially  in  areas 
of  high  stratification  or  curvature  in  the  main  tidal  channel 
[Dyer,  1974].  Irregularities  in  the  channel  edge,  side  creeks, 
and  intertidal  flats  can  also  contribute  to  along  channel 


5  of  16 


CO  1004 


MCKAY  AND  D1  IORIO:  VERTICAL  AND  HORIZONTAL  MIXING 


CO  1004 


dispersion  by  holding  water  inundated  on  flood  and  not 
releasing  it  back  into  the  main  channel  on  the  following  ebb, 
effectively  trapping  the  water,  and  returning  it  back  as  a 
different  water  mass  on  subsequent  ebb  tides  [ Fischer , 
1976]. 

[24]  As  shown  in  Figure  1 ,  our  deployment  in  the  lower 
Duplin  is  in  an  area  having  a  relatively  straight  channel  with 
few  side  creeks  or  irregularities.  The  width  of  the  channel  is 
approximately  B  =  200  m  with  the  eastern  side  approxi¬ 
mately  1  m  deeper  than  the  western  side  with  steep  sloping 
side  boundaries.  Following  the  work  of  Hansen  [1965]  and 
Fischer  [1976],  as  modified  by  Medeiros  and  Kjerfve 
[2005],  we  quantify  the  cross-sectional ly  integrated  and 
tidally  averaged  salt  flux  as 

d±l+'-j^aus)da  =  0,  (7) 

where  the  superscript  V  implies  a  volume  integration  of  the 
salinity  from  the  head  of  the  creek  where  the  salt  fluxes  are 
zero  to  the  measurement  location.  We  decompose  the  eross 
sectional  area  (a),  current  velocity  (w)  and  salinity  (s)  in 
terms  of  quantities  that  describe  the  residual  flow,  estuarine 
circulation  (£),  tidal  variability  (/),  vertical  shear  dispersion 
(K)  and  lateral  shear  dispersion  (T)% 

u(y,z%t)  =  U(z)  +  UE(z)  +  M,(0  +  Uy(zj)  +  UT(yj),  (8) 

a{t)  -  (H(t)  +  ht(i))B ,  (9) 

s(>\z,/)  =  S(t)  +  SE(z)  +  5,(0  +  Sy(z:t)  +  ST{y,t),  (10) 

where  the  cross  sectional  area  vanes  according  to  tidal  height 
and 


1  =  U  =  14,  UE  =  0 

(11) 

u,  =u  -  (2), 

(12) 

Uy  =  (u)  -  (2), 

(13) 

Ut  =  («)  -  (“)• 

(14) 

cation  is  much  greater  than  stratification  created  by  the 
estuarine  circulation  and  tidal  straining,  and  thus  we  can 
assume  that  the  lower  Duplin  during  our  measurement 
program  is  vertically  well  mixed  for  most  of  the  tidal  cycle. 
Substituting  equations  (8)  to  (10)  into  (7),  we  can  thus  write 
the  salt  flux  per  unit  width  of  the  estuary  as 

X-  +  Jj  J  [HUS  +  HUeSe  +  H(UySy) 

+  ~jj  f  h  (UTST)dy  +  U{u,s,)  +  U(h,s,) 

+  S{h,u,)  +  (h,u,s,)]dz  =  0.  (15) 

[26]  The  first  term  of  equation  (15)  is  the  storage  of  salt 
per  unit  width,  which  accounts  for  changes  in  salinity  and 
storage  volume.  The  first  term  within  the  integral  (HUS) 
represents  the  adveetion  of  salt  due  to  tidally  averaged 
residual  flows  reflecting  changes  in  freshwater  discharge, 
mean  depth  and  offshore  forcing.  The  next  term  is  the  net 
transport  of  salt  upstream  because  of  the  estuarine  circula¬ 
tion  that  moves  freshwater  down  stream  near  the  surface 
and  salt  water  upstream  at  depth.  The  next  two  terms  are  the 
vertical  and  transverse  shear  dispersion.  The  next  set  of 
terms  remaining  account  for  the  longitudinal  tidal  dispersive 
transport  due  to  correlations  between  tidally  varying  quan¬ 
tities  and  correspond  to  tidal  velocity-salinity  transport,  tidal 
salinity  transport,  tidal  wave  transport  (Stoke ’s  drift  which 
is  maximum  when  the  tide  resembles  a  progressive  wave 
and  is  zero  for  a  standing  wave)  and  the  tidal  pumping 
transport.  Rewriting  equation  (15)  as 

1  0(Sy)  7°  //  fh 

~g—gj - 1-  HUrS  =  -]  [UeSe  +  (UySy)]d:  —  —  J  {UrS-j)dv 

-  [U{u,s,)  +  Ur(h,s,)  +  S{h,u,)  +  (h,u,s,)\ 

(16) 

<i7> 

The  dispersive  flux  terms  can  then  be  parameterized  against 
the  along  channel  salinity  gradient,  dS/d: c,  to  obtain  an 
effective  along  channel  dispersion  coefficient,  Kx. 

3.2.1.  Estuarine  Circulation 

[27]  The  governing  equations  that  describe  the  residual 
flow  in  steady  state  are 


The  quantity  Ur  represents  the  mean  river  flow,? represents 
a  laterally  averaged  quantity  and  (ut)  =Uy  =  Ur  =  0.  Similar 
equations  can  be  written  for  the  salinity. 

[25]  As  the  instrumentation  was  located  in  the  center  of 
the  channel  and  the  bottom  salinity  measurement  was  fouled 
early  in  the  measurement  program  there  is  no  quantitative 
measurement  of  Uj{)\  /)>  and  s(y,  z).  As  a  result  we  cannot 
quantify  lateral  variations  in  the  longitudinal  fluxes  but  we 
expect  them  to  be  small  in  such  a  narrow  channel  with  steep 
sloping  sides  and  so  lateral  homogeneity  is  assumed.  In  order 
to  quantify  estuarine  circulation,  and  the  vertical  and  trans¬ 
verse  quantities  we  make  use  of  steady  state  formulae.  From 
the  energetics  of  mixing  (as  will  be  shown  below)  it  is 
expected  that  tidal  stirring  (mixing)  which  destroys  strati fi¬ 


0  = 


,  £ 

p  Ox  ‘  Oz2 


(18) 


OS  &Se 
UET<-K’-d?' 


(19) 


where  it  is  assumed  to  a  first-order  approximation  that  Az 
and  dpfdx  are  constant  with  depth.  The  barotropie  pressure 
gradient  associated  with  the  sea  surface  slope  (drj!dx)  drives 
the  fresh  water  out  on  top  and  the  baroelinic  pressure 
gradient  associated  with  the  along  channel  density  gradient 
(dpfdx)  drives  the  ocean  flow  in  at  depth. 
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[28]  Integrating  (18)  twiee  with  respect  to  z  and  applying 
the  boundary  conditions  that  there  is  no  stress  at  the  surface 
(A2dutdz\z=o  =  0)  and  no  flow  at  the  bottom  (w(z  =  -H)  =  0) 
gives 


To  find  the  sea  surfaee  slope  the  river  outflow  is  defined  as 
fL/fudz  =  HUr  and  as  such  the  depth-dependent  residual 
flow  is 


g  Op  II' 
p  dx  48/tr 


('-'©’■'©MK1'®! 


(21) 


-  UE(Z)  +  U(z)y 


(22) 


where  the  first  term  on  the  right  is  the  estuarine  velocity  and 
the  second  term  on  the  right  is  the  river  velocity. 

[29]  Substituting  the  estuarine  velocity  UE  into  (19), 
integrating  twice  with  respect  to  z,  applying  the  boundary 
conditions  of  no  flux  out  of  the  sea  surfaee  (A^S^dzl^  =  0) 
and  that  / fSfdz  =  0  gives  the  following  estuarine  salinity 
(for  K z  -  Az ): 

c  gflpdS  //5  (  1  I  (±\2  3/z\<  2/z\S\ 

E  pdxdxMA]\  12  2  \H)  4  \HJ  5  \n))‘ 

(23) 


The  depth  integrated  and  tidally  averaged  salt  flux  due  to 
vertical  shear  is  then 

f.  (27) 

whieh  is  also  always  negative  implying  a  net  salt  flux 
upstream. 

3.2.3.  Transverse  Shear 

[31]  In  a  similar  way  the  transverse  shear  dispersion  is 
calculated  for  a  channel  flow  of  u(y,  t)  -  U0(t)  (1  -  (ylb)2) 
where  B  =  2b  is  the  width  of  the  channel.  This  gives 
Ujiy,  t)  =  Un(t)(\l3  —  (ylb)2).  The  steady  state  equation  to 
solve  is 


(28) 


where  Ky  is  an  estimated  eross-estuary  mixing  parameter. 
Integrating  twiee  with  respeet  to  the  eross  channel  and 
applying  the  boundary  condition  of  no  fluxes  in  the  middle  of 
the  channel  (y  =  0)  because  of  symmetry  (KydSjJdy\}~o  =  0) 
and  that  ft  /,  Sjdy  =  0  gives 


Thus, 


f  j\uTST)dy=-{Ul) 


8  Hb2  OS 
945 ^  Ox  ’ 


(30) 


T  he  depth  integrated  salt  flux  due  to  the  estuarine  circulation 
is  then 


which  is  also  a  net  salt  flux  upstream.  The  difficulty  in  this 
formulation  is  quantifying  the  eross  channel  mixing  (Kv). 


UESEdz 


j9  /_g_ tips2 ir  os 

630  \  48p  dx )  A\  dx  ’ 


(24) 


which  is  relatively  constant  and  always  negative  implying 
that  there  is  a  net  salt  flux  upstream. 

3.2.2.  Vertical  Shear 

[30]  Bowden  [1965]  showed  how  the  interaction  between 
stratification  and  shear  in  the  mean  flow  can  give  rise  to 
horizontal  shear  dispersion.  If  the  velocity  profile  is  mod¬ 
eled  as  a  linear  function  with  depth  w(z,  /)  =  Ua(t)(\  +  (z/H)) 
then  Uv{z ,  /)  =  U0(t)(  1/2  +  z/H)y  where  the  surfaee  current  is 
defined  as  Ua.  The  steady  state  equation  to  solve  is  [Lewis, 
1 997] 


tPSy 


(25) 


which  represents  the  perturbed  state.  Integrating  with 
respeet  to  depth  twice  and  applying  the  boundary  conditions 
of  no  flux  out  of  the  surface  (KzdSyldz\^  =  0)  and  that 
J°  //vSy  dz  =  0  gives  the  salinity  (for  Kz  ~  Az) 


tll2dS(  I  l/z\2  l/z\A 
V  ~  °  Az  dx  (  24  "  4  (//)  4  6  (//)  )  ^  ^ 


3.3.  Energetics  of  Mixing 

[32]  The  important  processes  that  ean  lead  to  the  creation 
and  destruction  of  stratification  in  estuaries  are  from  the 
estuarine  circulation,  which  causes  fresh  water  outflow  at 
the  surface  and  oecanic  inflow  at  depth,  tidal  straining  of  the 
horizontal  density  field  by  the  vertical  current  shear,  whieh 
causes  an  increase  in  stratification  during  ebb  tide  and 
damps  the  turbulence  relative  to  the  flood  tide  [Simpson, 
1994;  Souza  and  Simpson ,  1997]  and  tidal  stirring  (mixing), 
which  is  a  percentage  of  the  turbulent  kinetie  energy  produc¬ 
tion  that  goes  directly  into  mixing  the  water  eolumn.  Stirring 
due  to  wind  and  waves  is  assumed  negligible  as  the  Duplin 
estuary  exists  on  the  back  side  of  a  eoastal  barrier  island. 

[33]  While  no  stratification  was  measured  at  the  mooring 
site  during  the  brief  period  when  both  surface  and  bottom 
measurements  were  available,  it  is  possible  that  some 
unmeasured  stratification  existed  at  later  times  and  this  eould 
have  an  impact  on  the  along  ehannel  salt  fluxes.  Unstratified 
estuaries  are  common  in  the  South  Atlantic  Bight  where 
relatively  high  tidal  energy  internets  with  relatively  low 
freshwater  fluxes  to  create  vertically  well  mixed  conditions 
[sce,c.g.,  Verity  etai,  1993;  Wiegertetal .,  1999].  Temporary 
stratification  ean  be  introduced  into  these  estuaries  through 
increases  in  freshwater  input  at  the  head,  through  an  increase 
in  river  or  groundwater  input  or  through  runoff  due  to  a  rain 
event,  as  well  as  the  decrease  in  tidal  energy  at  neap  tide. 
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While  these  stratification  events  are  likely  to  be  short-lived 
they  eould  have  an  impact  on  short-term  salt  fluxes. 

[34]  The  Duplin  River  has  no  regular  source  of  freshwater 
input;  its  freshwater  is  supplied  primarily  through  submarine 
groundwater  inputs  [Ragotzkie  and  Bryson,  1955]  and  local¬ 
ized  runoff  from  rain  events.  The  majority  of  the  ground- 
water  input  appears  to  enter  the  Duplin  upstream  of  the  lower 
tidal  excursion  [McKay,  2008],  and  as  a  result  the  along 
ehannel  salinity  gradient,  dS/dx  is  confined  to  the  lower 
Duplin  region.  In  this  seetion,  energetics  of  mixing  is  used  to 
see  if  there  are  times  when  stratification  can  be  created.  The 
conservation  of  density  and  its  depth  averaged  form  ean  be 
approximated  as 


dp  dp  d  — 


(31) 


d4 +«^=o, 

dt  <)x 


(32) 


assuming  no  fluxes  through  the  surface  and  bottom 
boundary  and  ffW  represents  the  turbulent  flux  of  mass  in 
the  vertical.  The  rate  of  change  of  potential  energy  anomaly 
is  then  defined  as  the  difference  between  the  mixed  and 
stratified  states  [Simpson  et  ai,  1990] 
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where  Rf  =  r/(T  +  1)  is  the  flux  Richardson  number 
representing  a  percentage  of  the  turbulent  kinetic  energy 
production  (calculated  as  P  =  —u'w'du/dz)  that  goes  into 
mixing  the  water  column,  creating  a  sink  of  turbulent 
kinetic  energy  (TICE)  defined  by  B  -  (g!p0)  ftvJ,  and  V  is 
the  mixing  efficiency.  The  first  term  on  the  right  hand  side 
of  (34)  represents  the  straining  that  ean  create  (positive  values 
on  ebb)  or  destroy  (negative  values  on  flood)  stratification; 
the  sceond  term  represents  the  destruction  of  stratification 
due  to  mixing  as  it  is  always  negative. 

4.  Results 

4.1.  Tidal  and  Residual  Properties 

[35]  Figure  3  (top)  shows  the  depth  distribution  of  the 
along  ehannel  velocity  u  (positive  on  ebb),  and  extrapolated 
to  the  surfaee  and  bottom  during  a  14  day  spring-ncap-spring 
transition  at  the  beginning  of  the  experiment.  Along  ehannel 
velocity  profiles  arc  generally  linear  with  a  minimum  near 
the  bed  and  a  maximum  near  the  surfaee.  Figure  3  (middle) 
shows  the  measured  near  surface  (5  meters  above  bottom 
(mab))  and  bottom  (1 .5  mab)  flow,  during  this  same  period. 
It  shows  more  clearly  the  spring/neap  cyelc,  the  ebb  domi¬ 
nance  of  flow  whieh  is  a  eommon  characteristic  of  this  class 
of  salt  marsh  estuary  [ Dronkers ,  1986;  Dyer ,  1997;  Blanton 
et  ai,  2002]  and  gives  an  indication  of  the  kind  of  shear  in 
the  water  eolumn.  This  shear  is  responsible  for  overturning 
the  water  column  and  keeping  it  essentially  well  mixed.  The 


entire  period  of  observation  is  shown  in  Figure  4,  separated 
into  residual  (40  hour  low-pass  filtered  shown  as  a  thiek  line) 
and  tidal  components  (shown  as  a  thin  line)  showing  the 
depth  averaged  quantities  for  the  along  ehannel  flow  with  the 
residual  flow  sealed  by  a  faetor  of  10  in  order  to  show  its 
variability  (Figure  4a);  showing  salinity,  whieh  is  taken  as 
representative  of  the  entire  water  column,  for  the  upper  and 
lower  Duplin  (Figure  4b);  and  showing  the  eenter  ehannel 
mean  and  tidally  varying  depth  (Figure  4e). 

[36]  A  strong  spring/neap  cycle  can  be  seen  in  both  along 
ehannel  tidal  velocity  and  tidal  depth.  Strong  spring  ebb 
flows  peak  as  high  as  1.2  m  s  1  while  weak  neap  ebb  flows 
are  as  slow  as  half  of  that,  and  the  tidal  height  ranges  from  1  m 
on  neap  to  2.6  m  on  spring  tide.  The  flow  is  ebb  dominant, 
showing  a  shorter,  stronger  ebb  which  contrasts  with  a  longer 
but  weaker  flood  This  is  due  to  frietional  distortion  of  the 
tide,  eaused  by  interactions  between  the  main  channel  and  the 
intertidal  areas,  whieh  ean  be  expressed  by  the  relationship 
between  the  M2  and  M4  tidal  constituents  [Dyer,  1997]. 
Residual  water  height  H  shows  some  inundation  that  may 
occur  when  the  along  and  eross  shelf  winds  are  strong. 

[37]  The  tidally  varying  salinity  varies  from  5  to  7  and  the 
mean  salinity  varies  between  23  and  28  in  the  lower  Duplin 
and  1 8  to  25  in  the  upper  Duplin.  The  decrease  in  salinity  on 
year  day  280  may  correspond  to  the  accumulated  precipi¬ 
tation  from  the  rain  events  shown  in  Figure  2  that  also 
causes  an  increased  net  outflow  This  pattern  is  also  seen  at 
the  peak  of  the  spring  tide  on  year  day  290  where  the 
salinity  decreases  when  there  is  a  net  outflow. 

[38]  The  tidally  averaged  current  flow  is  shown  in  Figure  5 
along  w  ith  the  theoretical  estuarine  circulation  defined  in  (21 ) 
(fora  value  of  Az-  1.1  x  10”3  m"/s,<9p/<9x  =  2.1  x  10~4kgW 
and  H  =  6.5  m).  The  resulting  river  flow  Ur  =  .021  m/s  gives 
the  best  fit  to  the  measured  profile. 

4.2.  Reynolds  Stresses  and  TKE  Production 

[39]  In  Figure  3  (bottom)  we  show  the  stress  ( rjp  =  —  u'W) 
through  the  water  eolumn  during  the  first  spring-ncap-spring 
period.  Maximum  stresses  originate  near  the  bed  and  prop¬ 
agate  upward  into  the  water  column  decreasing  in  the  upper 
layers.  The  periods  of  highest  stress  are  found  during  the 
spring  tide  centered  around  maximum  flood  and  ebb  when 
the  stresses  reach  as  high  as  ±1 .2  x  10  3  m2  s  near  the  bed 
and  show  an  M4  distribution.  There  is  a  pronounced  tidal 
asymmetry  in  the  stress  production.  They  are  lower  on  flood 
than  ebb  and  the  stresses  on  ebb  penetrate  significantly 
further  into  the  water  column  than  on  flood,  due  to  peak 
ebb  velocities  being  much  higher,  though  of  shorter  duration, 
than  flood.  At  neap  tide  there  is  greatly  decreased  turbulent 
stress  on  maximum  ebb  and  flood  with  consequently  less 
penetration  of  the  stress  into  the  water  eolumn.  Near  slaek 
water,  at  both  spring  and  neap  tides,  the  stresses  are  uniformly 
low  and  patehy  throughout  the  water  column. 

[40]  Figure  4d  shows  a  time  series  of  rjp,  centered 
1.5  mab,  for  the  entire  deployment  period.  Both  fortnightly 
and  semidiurnal  variability  in  stress  production  arc  elearly 
visible  with  peak  stresses  during  spring  tides  being  much 
higher  than  during  neap  tides,  and  with  the  peak  stresses 
during  the  second  spring  tide  being  more  than  double  those 
during  the  first.  Semidiurnal  variability  is  seen  in  that  peak 
positive,  ebb  tide,  stresses  are  significantly  higher  than  peak 
negative,  flood  tide,  stresses  due  to  the  asymmetry  in  ebb  and 
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Figure  4.  Tidally  averaged  (40  hour  low-pass  filtered  shown  as  thick  lines)  and  tidally  varying  (shown 
as  thin  lines)  quantities  for  (a)  depth  averaged  along  channel  velocity  (U  is  sealed  by  a  factor  of  10), 
(b)  surface  salinity  for  the  lower  (solid  line)  and  upper  (dashed  line)  Duplin,  and  (e)  center  channel  water 
depth,  (d)  The  Reynolds  stress  rjp  and  (e)  the  TKE  production  (P)  both  at  1.5  mab  are  shown  together 
with  spnng  (S)  and  neap  (N)  tide  labels. 


flood  velocities.  Cross  channel  stresses,  rjp  —  —  v'w',  are 
not  shown  in  these  plots  as  they  are  uniformly  low,  owing 
to  the  low  cross  channel  velocity. 

[41]  Figure  4e  also  shows  the  turbulent  kinetic  energy 
production  ( P )  at  1 .5  mab.  The  depth  variability  (not  shown) 
mirrors  the  stress  distribution  as  the  velocity  profile  is  very 
close  to  linear  (as  will  be  shown).  Production  of  TKE  is 
maximum  on  spnng  tide  and  maximum  at  the  bottom 
decreasing  upward  into  the  water  column  Peak  TKE  pro¬ 
duction  is  approximately  5  x  10  4  W  kg”1. 

4.3.  Vertical  Mixing 

[42]  From  (2)  the  vertical  eddy  difTusivity  of  momentum, 
AZi  is  calculated  using  measurements  of  the  turbulent  stress, 
rrjp  =  — w'w'  and  vertical  shear,  duldz.  Sample  velocity 
profiles  from  maximum  flood  to  maximum  ebb  on  18 
October  (year  day  291)  at  the  peak  of  the  strong  second 
spnng  tide  are  shown  in  Figure  6.  The  velocity  profiles  can  be 
seen  to  be  nearly  linear  through  the  portion  of  the  water 
column  resolved  by  the  ADCP  with  little  curvature,  thus 
having  a  nearly  constant  shear.  The  log  layer  fit  below  the 


lowest  ADCP  bin  is  also  shown.  As  the  shear  is  nearly 
constant  with  depth  the  profile  of  Az  over  our  measurement 
region  follows  a  near  linear  distnbution  with  depth  What 
parabolic  shape  there  is  to  the  Az  profile  would  occur  below 
the  lowest  measured  ADCP  bin  where  the  log  layer  pre¬ 
sumably  applies.  While  stresses  arc  lower  on  flood  than  on 
ebb,  shear  is  less  on  flood  than  on  ebb 

[43]  Significant  stresses  are  mainly  generated  around  the 
time  of  maximum  flood  and  ebb  at  spring  tide  and  these  will 
thus  be  the  times  of  the  greatest  vertical  mixing.  The  depth 
normalized  profiles  of  Az  at  the  times  of  maximum  flood 
and  ebb  at  the  peak  of  the  spring  tide  is  shown  in  Figure  6. 
Like  the  stresses,  the  values  are  greatest  near  the  bed,  with 
values  of  1.2  x  10“2  m2  s_1  rapidly  decreasing  with  height 
above  bottom  to  near  zero  close  to  the  surface.  Flood  and 
ebb  values  are  very  similar  and  both  follow  the  plotted  least 
squares  linear  fit,  which  has  an  r2  =  0.81. 

[44]  Neap  tide  stresses  are  similarly  concentrated  around 
times  of  maximum  flood  and  ebb  but,  as  can  be  seen  in 
Figure  3,  the  periods  of  significant  stress  generation  are 
mueh  shorter  than  on  spring  tide  and  the  high  stress  region 
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Figure  5.  The  estuarine  circulation  and  river  flow  as 
modeled  by  equation  (21)  compared  to  the  measured  tidally 
averaged  flow. 

does  not  penetrate  nearly  as  far  up  into  the  water  column. 
When  the  stresses  are  low,  patchy  and  variable,  they  do  not 
present  a  well  defined  profile  in  Az. 

[45]  As  water  masses  mix  through  the  water  column  they 
experience  the  entire  vertical  range  of  mixing  levels  and 
thus  a  depth  averaged  mixing  quantity  is  approximated  by 
Kz  ~  Az  which  is  shown  in  Figure  7  (top).  Both  tidal  and 
fortnightly  vanability  arc  apparent  with  values  fluctuating 
strongly  during  the  tidal  cycle  tied  to  the  generation  of  stress 
at  the  bed.  Az  shows  M4  vanability  tied  to  the  two  periods 
of  maximum  flood  and  ebb  during  each  lunar  day.  The 
40  h  low-pass  filtered  value  of  Azy  shown  by  the  thick  line, 
varies  on  a  fortnightly  time  scale  being  highest  on  spring 
tides,  centered  around  year  day  277  and  290,  in  the  range  of 
4  6x10  1  m2  s  1  and  lowest  on  neap  tide,  being  half  to  a 
third  of  the  spring  tide  values.  Even  the  low  values,  however, 
represent  significant  vertical  mixing  in  this  shallow  system 
[Lewis,  1997] 

[46]  Equation  (6)  is  used  to  calculate  the  approximate 
time  for  complete  vertical  mixing  during  a  tidal  period, 
using  the  measured  mean  water  depth  (//),  and  the  tidally 
and  depth -averaged  value  of  Az.  This  is  shown  in  Figure  7 
(bottom)  plotted  in  minutes.  Even  at  its  slowest,  complete 
vertical  mixing  is  achieved  in  slightly  over  1  hour  and  at 
spring  tide  takes  as  little  as  20  minutes.  This  is  very  much 
less  than  the  tidal  period  and  implies  that  the  lower  Duplin 
stays  well  mixed  at  all  times  unless  acted  upon  stratifying 
influences  like  tidal  straining  or  the  estuarine  circulation. 

4.4.  Stratifying  Influences 

[47]  While  the  rain  events  and  groundwater  input  could 
serve  to  create  stratification  through  straining  processes, 
thus  potentially  modifying  the  along  channel  salt  fluxes,  the 
tidal  mixing  (stirring)  serves  to  homogenize  the  water 
column.  Figure  8  shows  the  energetics  of  stirring  and 


straining  according  to  equation  (34)  in  order  to  determine  if 
there  were  times  when  stratification  would  become  impor¬ 
tant.  The  depth  averaged  turbulent  kinetic  energy  ( P )  was 
calculated  and  a  flux  Richardson  number  of  Rf  =  0.05  is 
assumed  Di  Iorio  and  Kang  [2007]  measure  the  mixing 
efficiency  T  in  the  Altamaha  River  estuary  which  is  partially 
to  well  mixed  and  show  that  it  is  on  average  0.05  but  can 
approach  the  upper  limit  of  0.2.  As  expected  significant 
mixing  occurs  during  spring  tide  which  is  greater  than  the 
stratifying  influences  of  straining. 

[48]  The  straining  term  is  calculated  from  the  measured 
velocity  profiles  and  the  along  channel  density  gradient  is 
calculated  from  the  tidally  averaged  along  channel  advcc- 
tion  of  mass  past  the  mooring  as  dpldx  =  (—  1  !U)dp!dt.  The 
straining  effect  shows  that  stratification  is  created  on  ebb 
tide  (positive  values)  and  destrati fieation  on  flood  (negative 
values)  with  greater  effects  on  ebb  owing  to  its  greater  flow 
and  shear.  Combining  the  stirring  and  straining  effect,  the 
rate  of  change  of  the  potential  energy  anomaly  is  always 
negative  except  for  brief  periods  during  neap  tide  when 
stratifying  influences  come  into  play.  These  periods  of  strat¬ 
ification  are  very  short  and  any  stratification  that  exists  at  the 
slack  tide  will  be  quickly  eroded  away  during  the  subsequent 
flood.  1  f  /?/  ~  0.008  then  there  can  be  stratifying  tendencies 
on  every  ebb  and  if/?^>  0.05  then  no  stratification  will  exist. 
According  to  Geyer  et  al.  [2008]  a  distribution  of  /^values 
were  obtained  in  the  Merrimack  estuary  during  an  ebb  tide 
with  a  median  value  of  0. 15  and  with  low  values  <0.01  during 
times  of  weak  stratification.  As  a  result  the  Duplin  River  can 
be  assumed  to  be  essentially  well  mixed  or  weakly  stratified 
only  during  neap  tides, 

4.5.  Horizontal  Transport 

[49]  The  measured  salt  fluxes  per  unit  width  from  the 
various  tidal  correlations  in  (16)  arc  shown  in  Figure  9a  as 
well  as  their  net  sum  (thick  solid  line).  The  largest  of  these 
tidal  correlation  terms  is  the  salt  flux  into  the  Duplin  due  to 
the  correlation  between  tidal  velocity  and  depth  (dashed  line). 
The  deeper  water  on  the  weaker  but  longer  flood  allows  more 
salt  to  flow  in  than  the  shallow  waters  on  the  stronger  but 
shorter  ebb  allow  out.  The  next  largest  term  is  a  net  flux  out 
due  to  the  correlation  between  velocity  and  salinity  (dotted 
line)  caused  by  the  higher  ebb  tide  velocity  exporting  salt 
which  has  been  mixed  in  on  flood.  Both  of  these  terms  show  a 
strong  spring/neap  modulation  with  each  approaching  zero 
on  neap  tide  and  being  greatest  on  spring.  Of  the  remaining 
terms  the  correlation  between  depth  and  salinity  (thin  line)  is 
negated  by  being  multiplied  by  the  very  low  mean  velocity, 
on  the  order  of  ±1  cm  s  \  and  the  triple  correlation  is 
similarly  small.  The  net  tidal  flux  is  shown  by  the  thick  line 
as  a  net  flux  of  salt  into  the  Duplin,  where  it  follows  the 
same  strong  spring/neap  modulation  as  do  its  components. 

[50]  The  estimated  fluxes  due  to  the  estuarine  circulation, 
vertical  and  transverse  shear  dispersion  calculated  from 
equations  (24),  (27)  and  (30),  respectively,  are  shown  in 
Figure  9b.  Of  these  three  effects  the  vertical  shear  dispersion 
is  the  most  dominant  with  values  approaching  —0.2  m2/s  and 
longitudinal  mixing  coefficients  KxV  ~  60  m/s.  It  should  be 
noted  that  the  linear  velocity  profile  used  is  actually  an 
overestimation  of  the  kind  of  shear  that  actually  persists. 
The  effect  of  the  estuarine  circulation  is  a  factor  of  104 
smaller.  For Ky=  [0. 1  0.5  10  1 0]  m2/s,  cross  channel  salinity 
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t=0h  t=+1h  t=+2h  t=+3h  t=+4h  t=+5h  t=+6h 


Figure  6.  (top)  Profiles  of  along-channel  velocity  through  one  tidal  cycle  from  maximum  ebb  to 
maximum  flood  on  year  day  291.  (bottom)  Profiles  of  Az  at  maximum  ebb  and  maximum  flood,  during 
the  times  of  significant  stress  generation  on  spring  tide.  Depths  are  normalized  with  respect  to  channel 
depth.  The  linear  fit  in  the  Figure  6  (bottom)  has  an  r2  of  0.81. 


Year  Day  2005  (275  =  2  OCT) 

Figure  7.  Time  series  of  depth  averaged  vertical  eddy  diffusivity  of  momentum  (A2)  and  vertical  mixing 
time  (/)  during  the  deployment  period.  The  thick  lines  indicate  40  hour  low-pass  tidally  averaged  values. 
The  labels  S  and  N  indicate  times  of  spring  and  neap  tide,  respectively. 
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Figure  8.  Energetics  of  (top)  mixing  and  (middle)  straining  for  the  destruction  (negative  values)  and 
creation  (positive  values)  of  stratification  together  with  (bottom)  their  sum. 


gradients  correspond  to  a  change  in  salinity  of  [4.17  0.83 
0.42  0.04]  in  100  m,  respectively.  It  is  unlikely  that  a  cross 
channel  salinity  gradient  of  4.2  in  100  m  exists  in  this 
environment  and  so  we  limit  the  value  of  Ky  >  0.1  m2/s 
and  assume  that  Kv  ~  Kxy.  The  salt  flux  from  transverse 
shear  dispersion  calculated  in  this  way  is  a  factor  of  100 
smaller  than  the  vertical  shear  dispersion.  If  the  transverse 
and  vertical  shear  dispersions  are  to  be  of  the  same  order 
then  Ky  ^  1.0  m2/s. 

[51]  As  all  these  estimated  salt  fluxes  are  very  much 
lower  than  tidal  dispersion  terms,  and  thus  the  effects  of 
possible  stratification  and  hence  shear  dispersion  are  judged 
to  be  minor  players  in  the  long-term  salt  budget  during  this 
time.  Uncles  et  al.  [1985]  also  showed  that  in  an  estuary  of 
the  same  scale  as  the  Duplin  that  the  transverse  shear  dis¬ 
persion  was  much  smaller  than  the  vertical  shear  dispersion. 
The  dominant  driver  of  along  channel  horizontal  dispersion  is 
then  not  shear  acting  on  salinity  gradients,  but  the  tidal 
correlations,  which  act  much  more  slowly  but  with  greater 
energy  and  on  larger  scales.  Jay  and  Smith  [  1 990]  outline  that 
barotropie  adveetion  of  salt  dominates  over  density-driven 
circulation  in  maintaining  the  salt  balance  in  weakly  stratified 
estuaries. 

[52]  Figure  9e  shows  the  net  tidal  flux  term  (dashed  line) 
along  with  the  flux  of  salt  due  to  residual  adveetion  and  the 
change  in  mean  water  depth  tied  to  changes  in  storage  volume 
(thin  solid  line).  Despite  the  low  residual  velocity  the  residual 
adveetive  term  is  much  greater  than  the  tidal  flux  term  due  to 
the  high  value  of  the  mean  salinity  compared  to  its  tidal 
fluctuations.  Net  adveetion  varies  strongly  on  a  spring/neap 


cycle  showing  a  pulse-like  export  of  salt  on  spring  tides 
followed  by  a  weaker,  but  longer,  period  of  import  of  salt  on 
neap  tides.  This  pattern,  which  has  been  observed  before  in 
well  mixed  to  partially  stratified  estuaries  [sec,  e.g.,  Nunes 
and  Lennon,  1987;  Bow  en  and  Geyer ,  2003]  is  believed  here 
to  be  due  not  only  to  tidal  modulation  of  upstream  dispersive 
mechanisms  but  also  to  variations  in  the  flux  of  fresh 
groundwater  into  the  upper/middle  Duplin  caused  by 
changes  in  the  tidal  range  on  the  spring/neap  cycle  [Schultz 
and  Ruppel ,  2002]. 

[53]  On  neap  tide,  when  the  tidal  range  is  small,  we 
hypothesize  that  little  fresh  groundwater  is  pumped  into 
the  uppcr/middle  Duplin  and  thus  little  pressure  is  exerted 
on  water  in  the  lower  Duplin  and  residual  adveetion  nears 
zero.  Tidal  fluxes  then  bring  salt  into  the  lower  Duplin  and 
the  salinity  in  this  region  increases.  With  the  increased  tidal 
range  at  spring  tide,  more  fresh  groundwater  is  pumped  into 
the  upper/middle  Duplin  which  then  flows  down  to  the 
lower  Duplin.  Residual  adveetion  increases  as  this  freshwater 
mass  forces  saltier  water  out  of  the  lower  Duplin  and  salinity 
decreases  as  it  mixes  into  the  lower  Duplin  water  mass.  Tidal 
dispersion  of  salt  into  the  lower  Duplin  increases  with  the 
increased  tidal  energy  but  it  is  not  enough  to  overcome  the 
diluting  effects  of  the  fresh  groundwater.  When  compared 
to  lower  Duplin  salinity  in  Figure  4b  the  net  flux  can  be 
seen  to  closely  correlate  with  changes  in  the  mean  salinity 
in  the  lower  Duplin.  The  sharp,  but  short-term,  decrease  in 
salinity  around  year  day  280  as  well  as  the  lesser,  but  longer 
lasting,  decrease  around  year  day  290  both  correlate  with 
spring  tide  induced  net  export  as  shown  in  Figure  9d  while 
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Year  Day  2005  (275  =  2  OCT) 


Figure  9.  (a)  Tidal  pumping  salt  fluxes  together  with  their  sum.  (b)  Estimated  salt  fluxes  due  to  the 
estuarine  circulation  and  due  to  vertical  and  transverse  shear  dispersion,  (e)  Adveeti ve  and  dispersive  fluxes 
and  their  sum  showing  the  net  salt  flux  per  unit  width,  (d)  A  comparison  between  the  storage  of  salt  with  the 
net  salt  flux. 


the  slow'  recovery  periods  correlate  with  the  low  net  influx  of 
salt  on  neap  tide. 

[54]  Figure  9d  also  shows  the  estimated  longitudinally 
and  depth  integrated,  and  tidally  averaged  storage  of  salt, 
per  unit  w  idth,  along  with  the  negative  flux  of  salt,  per  unit 
width,  past  the  mooring.  While  it  can  be  difficult  to  accurately 
estimate  upstream  salt  storage  [Simpson  et  ai,  2001],  a 
general  estimation  may  be  made  by  combining  salinity 
measurements  with  an  estimate  of  the  time  varying  volume 
of  water  upstream  as  estimated  from  tidal  height  and 
bathymetry.  This  estimate  will  not  be  precise  but  will  show 
the  general  magnitude  and  trend  of  the  upstream  storage  of 
salt. 

[55]  Both  the  sign  and  the  shape  of  the  curves  compare 
reasonably  well.  The  general  trend  of  the  curves  suggest 
that  salinity  generally  increases  during  neap  tide  followed 
by  a  decrease  in  salinity  during  spring  tide.  This  is  believed 
due  to  fresh  groundwater  transported  down  from  the  middle 
Duplin  area. 

[56]  To  get  an  along  channel  dispersion  coefficient  Kx 
for  the  lower  Duplin  we  normalize  the  dispersive  fluxes 


(dominated  by  tidal  processes)  by  the  mean  water  depth  and 
the  along  channel  residual  salinity  gradient,  dS/d; c,  shown  in 
Figure  10.  Figure  10  (top)  shows  the  depth  normalized  tidal 
dispersion  term  from  Figure  9  while  the  Figure  10  (middle) 
shows  the  along  channel  residual  salinity  gradient,  calculated 
from  the  tidally  averaged  along  channel  advection  of  salt 
past  the  mooring  as  dS/dx  =  (~\/U)dS/dt.  The  salinity 
gradient  is  positive,  showing  saltier  water  toward  the  mouth 
of  the  Duplin,  with  its  greatest  value  on  neap  tide  as  decreased 
tidal  energy  pumps  less  salt  into  the  upper  reaches  of  the 
creek.  The  ratio  of  the  depth  normalized  tidal  pumping  flux 
with  the  salinity  gradient  gives  a  value  for  the  along  channel 
tidal  dispersion  coefficient  Kx.  This  is  shown  in  Figure  10 
(bottom)  and  varies  strongly  on  a  fortnightly  cycle  from 
near  1 000  m2  s  1  on  spring  tide  to  close  to  zero  on  neap  tide 
with  a  long-term  mean  of  approximately  500  m2  s  *. 

5.  Discussion  and  Conclusions 

[57]  We  have  presented  measurements  near  the  mouth  of 
the  Duplin  River,  which  is  a  small  and  generally  vertically 
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Figure  10.  Time  series  of  the  dispersive  salt  fluxes  per  unit  width  normalized  by  mean  depth,  the  along 
channel  salinity  gradient,  and  the  resultant  along-channel  dispersive  coefficient  (Kx)  during  the  deployment 
period 


well-mixed  tidal  creek.  This  creek  system  serves  as  the 
main  conduit  between  Doboy  Sound  and  the  intertidal 
marshes  of  the  Sapclo  Island  National  Estuarine  Research 
Reserve.  The  measurements  described  here  directly  resolve 
the  vertical  velocity  structure  through  most  of  the  6.5  m 
mean  water  depth  of  the  lower  Duplin  as  well  as  allow  the 
estimation,  through  the  variance  method,  of  the  stresses 
which  give  rise  to  vertical  mixing.  Through  salt  budget 
methods  we  also  quantify  the  along  channel  tidal  dispersion. 
Spanning  a  34  day  period,  consisting  of  two  spring  and  two 
neap  tides  along  with  the  transitions  between  these  tides,  our 
measurements  resolve  the  wide  range  of  time  scales,  from 
turbulent,  to  tidal  and  fortnightly,  which  influence  the  mixing 
in  these  waters. 

[58]  Estimates  of  the  vertical  turbulent  stress,  rjp,  show 
modulation  on  both  tidal  and  fortnightly  frequencies  with 
maximum  stress  being  generated  near  the  bed  on  maximum 
flood  and  ebb  at  spring  tide  and  linearly  decreasing  with 
height  above  bottom.  Neap  tide  stresses  are  weaker,  concen¬ 
trated  near  the  bottom  and  have  no  variability  with  depth. 
Internally  generated  turbulence  plays  a  larger  role  near  times 
of  slack  water  with  stresses  being  more  broadly  distributed 
through  the  water  column. 

[59]  Due  to  the  nearly  linear  profile  of  velocity  with  depth 
and  the  constant  density  the  Richardson  number  is  almost 
always  much  less  than  the  stability  criteria  of  0.25  and  thus 
from  the  diffusion  of  momentum  we  can  approximate  the 
vertical  mixing  of  salt  using  a  turbulent  Prandtl  scaling  of 
~1 .  From  this,  we  estimate  that  the  time  for  complete 


vertical  mixing  vanes  from  20  30  minutes  on  spring  tide 
to  40  60  minutes  on  neap  tide  and  is  at  all  time  very  much 
less  than  the  tidal  period,  thus  ensunng  that  the  water  column 
stays  well  mixed  at  all  times.  A  comparison  of  energetics  for 
mixing  (Rf  =  .05)  with  strained  induced  periodic  stratification 
show  shows  that  only  for  brief  periods  dunng  neap  tide  can 
stratification  exist. 

[60]  Tidal  dispersive  fluxes  are,  by  nature,  tidally  averaged 
so  we  cannot  comment  on  intertidal  variability  and  time 
scales.  These  fluxes  show  a  strong  spring/ncap  cycle  with 
maximum  upstream  dispersive  fluxes  being  found  at  spring 
tide  while  neap  tide  fluxes  approach  zero  but  still  being  much 
greater  than  theoretically  predicted  vertical  and  transverse 
shear  effects  and  effects  due  to  the  estuarine  circulation.  The 
upstream  flux  due  to  the  correlation  between  the  tidally 
varying  depth  and  velocity  (tidal  wave  transport)  is  the 
primary  driver  of  along  channel  dispersion,  but  it  is  some¬ 
what  counterbalanced  by  an  opposing  dispersive  flux  due 
to  the  correlation  between  tidally  varying  velocity  and 
salinity.  These  results  are  similar  to  those  obtained  from 
Kaluarachchi  et  al.  [2003]  with  the  Stokes  term  being  the 
most  significant  upstream  transport  of  salt  in  Newark  Bay 
and  Arthur  Kill  transects  but  since  their  model  contains  only 
one  grid  across  the  channels  they  could  not  quantify  trans¬ 
verse  shear  effects.  These  dispersive  salt  fluxes  are  generally 
overwhelmed  by  the  net  advective  flux  of  salt  due  to  the 
residual  flow  out  of  the  creek.  Residual  flows  in  this  envi¬ 
ronment  can  be  tied  to  changes  in  groundwater  input,  and 
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changes  in  rainfall  and  oceanic  forcing  from  meteorological 
events. 

[61 ]  Interaction  between  the  main  along  channel  flow  and 
the  channel  edge  may  result  in  enhanced  along  channel 
dispersion  due  to  the  trapping  of  water  in  dead  zones,  side 
creeks  and  embayments,  and  intertidal  marshes  and  mudflats 
[Fischer  et  ai ,  1979]  This  water  is  diverted  from  the  main 
flow  on  flood  and  rejoins  it  on  a  different  stage  of  ebb,  or 
on  a  different  ebb  altogether,  where  it  mixes  with  different 
water  masses  both  creating  cross  channel  gradients  and 
enhancing  (or  retarding)  effective  along  channel  dispersion. 
Trapping  in  dead  zones,  generally  eddies  associated  with 
channel  bends  or  small  embayments  and  obstructions  at  the 
channel  edge,  is  one  of  the  more  well  studied  contributors 
to  enhanced  dispersion  due  to  trapping.  A  number  of  theo¬ 
retical  and  numerical  studies,  generally  in  nontidal  rivers, 
[sec,  e.g.,  Pumama ,  1988;  Davis  et  ai ,  2000;  Davis  and 
Atkinson ,  2000;  Lees  et  ai,  2000]  have  established  that  in 
areas  of  high  curvature  or  channel  edge  irregularities  this  can 
be  a  major  factor  in  enhancing  along  channel  dispersion, 
though  the  exact  amount  is  highly  dependent  on  the  local 
morphology. 

[62]  As  can  be  seen  in  Figure  1,  our  deployment  is  in  an 
area  with  a  straight  channel  with  few  side  creeks  or 
irregularities,  which  results  in  few  dead  zones  to  act  on 
along  channel  dispersion.  The  nearest  significant  side  creek, 
Bam  Creek,  is  located  nearly  an  entire  tidal  excursion  from 
the  mooring  location,  thus  making  it  unlikely  that  it  plays  a 
role  in  enhancing  dispersion  through  tidal  trapping  in  the 
lower  Duplin.  While  there  arc  intertidal  marshes  bordering 
the  mooring  site  there  is  also  extensive  hard  upland, 
associated  with  Sapclo  Island  and  several  marsh  hammocks, 
or  isolated  areas  of  hard  upland,  which  have  been  shown  to 
greatly  reduce  the  effects  of  the  intertidal  areas  on  local 
processes,  especially  as  compared  to  areas  further  upstream 
[McKay  and  Di  Iorio,  2008].  Due  to  these  considerations  it 
is  reasonable  to  neglect  the  contributions  of  channel  edge 
effects  and  tidal  trapping  to  the  along  channel  dispersion  for 
the  lower  Duplin. 

[6i]  The  net  water  flux  is  generally  out  at  spring  tide.  The 
increased  tidal  range  on  spring  tide  presumably  increases 
the  hydraulic  head  acting  on  the  aquifer  increasing  the 
groundwater  flux  into  the  middle  and  upper  Duplin  [Schultz 
and  Ruppel,  2002].  This  increases  the  downstream  pressure 
gradient  and  thus  increases  the  net  outflow  of  water.  On  neap 
tide,  groundwater  input  into  the  upper  and  middle  Duplin 
decreases  and  tidal  pumping,  which  is  no  longer  opposed  as 
strongly  by  residual  export,  brings  salt  in  from  Doboy  Sound. 

[64]  The  total  salt  flux  due  to  tidal  dispersive  and  advcc- 
tivc  fluxes  correlates  well  with  the  trends  in  lower  Duplin 
salinity  during  the  deployment  period.  Examining  the  net 
horizontal  flux,  it  becomes  apparent  that  the  Duplin  shows  a 
net  export  of  salt  for  a  brief  period  on  spring  tide,  which  is 
followed  by  a  longer  penod  of  cither  net  import  or  near 
neutral  conditions.  Salt  export  is  then  in  a  pulsating  pattern 
tied  to  the  fortnightly  spring/ncap  tidal  frequency. 
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